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Abstract. - We study a mesoscopic superconducting ring threaded by a magnetic flux when 
the single particle level spacing is not negligible. It is shown that, for a superconducting ring 
with even parity, the behavior of persistent current is equivalent to what is expected in a 
bulk superconducting ring. On the other hand, we find that a ring with odd parity shows 
anomalous behavior such as fluxoid quantization at half-integral multiples of the flux quantum 
and paramagnetic response at low temperature. We also discuss how the parity effect in the 
persistent current disappears as the temperature is raised or as the size of the ring increases. 



What happens to superconductivity when the sample is made very small? Anderson ffl already 
addressed this fundamental question in 1959 and argued that as the size of a superconductor 
decreases and, accordingly, the average level spacing 5 becomes larger than the BCS gap A, 
superconductivity is no longer possible. Recent experiments on ultrasmall "superconducting" 
nanoparticles [0 have led to reconsider this old, but fundamental question. In a series 
of experiments the authors of Q studied transport through nanometer-scale Al grains and 
succeeded to get the discrete eigenspectrum of a single superconducting grain. The results 
were found to depend on the parity, i.e. on the electron number in the grain being even or 
odd. These experiments initiated several theoretical investigations, von Delft et al. [|| used a 
model of uniform level spacing in a parity-projected mean field theory Q| and found that the 
breakdown of superconductivity occurs at a value of <5/A ~ 0(1) which is parity-dependent. 
This parity effect has been shown to increase when including the effects of level statistics || . 
Effects of quantum fluctuations ||, canonical description of BCS superconductivity 0, as well 
as transport theory for a nanoparticle coupled to superconducting leads Q have also been 
subjects of the study along this line. 

Another interesting example for studying the size effect on superconductivity is a mesoscopic 
superconducting ring threaded by a magnetic flux $. It is well known that a conventional BCS 
superconducting ring exhibits fluxoid quantization at integer multiples of the flux quantum 
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$o = hc/2e and a diamagnetic response at <£> = n$o with n being an integer ||. In this Letter, 
wc address the following question which is essentially the same as in a simply connected grain: 
What happens in a superconducting ring when the size of the ring becomes very small? For 
this purpose, we adopt the parity-projected mean field theory [Q for an ideal mesoscopic ring. 
It is shown that the order parameter strongly depends on the parity when the size of the 
ring is small enough, as in the case of a grain. The most dramatic feature we show in our 
study is the behavior of the supercurrent (or persistent current) which strongly depends on the 
parity. For a ring with even parity, the behavior of the supercurrent is identical to that of a bulk 
superconducting ring. It exhibits fluxoid quantization at integer multiples of the flux quantum, 
and a diamagnetic response for small deviations of the flux from the integer multiples of the 
flux quantum. On the other hand, the characteristics are found to be very different for a small, 
odd parity ring. In a superconductor with odd parity, there is one unpaired quasiparticle. We 
find that at low temperature the existence of this quasiparticle drives the superconducting ring 
to a half-integral fluxoid quantization and a paramagnetic response at small values of the flux. 
We further show that the anomalous behavior in an odd-parity superconductor disappears at 
temperatures higher than the level spacing of single electron spectra, where one recovers the 
behavior of a conventional superconducting ring. Finally, this parity-dependent behavior of 
supercurrent is shown to disappear in the thermodynamic limit. 

The ideal superconducting ring can be described by the Hamiltonian: 

j<r i,j 

The single particle energy is given by e° = 2 ^ fl2 (J — f /2) 2 , where R is the radius of the ring, 
f = $/<i>Q is the external flux divided by the flux quantum $o = hc/2e, and j is an integer 
which corresponds to an angular momentum quantum number. This is obtained by solving the 
Schrodinger equation in a ID noninteracting ring. Note that j = —j. A is the dimcnsionless 
BCS coupling constant. 5 — H 2 N/8mR 2 is the level spacing at the Fermi energy, where N is 
the number of electrons in the ring. We don't take into account the Zeeman splitting, namely 
h, because it is negligible unless the radius of the ring is very small. For $ ~ $o with a uniform 
magnetic field, the ratio of the level spacing 5 to the Zeeman splitting is proportional to R 
which is estimated as 5/h ~ l0 20 r s R/m 2 , where r s is the average distance between electrons. 
For example, in a typical superconductor such as Al with R ~ 1/im, 5/h ~ 10 4 . 

A simple way of describing a mesoscopic superconductor with fixed number parity P 
(denoted by e for even, and o for odd parity) is to adopt the parity-projected grand canonical 
partition function Q 

Z P {n)=Tth L l±{-l) N ]e-K H -» N \ (2) 

We evaluate Zp using the BCS-type mean field approximations, which consists in neglecting 
quadratic terms of the fluctuations : 

- ( c h4i) c n c n + 4i4i( c ji c n) - (4Ai^ c h c n) ( 3 ) 

+ % (4 t c *t (4i°i0 + <4r c *t) c li c a - (4^^)(4i c n)) ■ 

The ensemble average (• ■ •) should be evaluated in a given parity P — e or P = o. The 
first three terms on the r.h.s of Eq.([|) correspond to the mean-field approximation for the 
superconducting pairing. The last three terms are usually not considered since they give 
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no contribution in the thermodynamic limit. However, those terms cannot be ignored in 
mesoscopic systems, though they were neglected in the previous mean-field description for 
ultrasmall superconducting grains Note that the validity of our mean- field treatment is 
limited to the 5 < A limit where the usual BCS approximation can be applied. As a result we 
get the following expression for the mean-field Hamiltonian 

H = C P + J2Ej7llja+^N , (4) 

3° 

where jj a (jj a ) destroys (creates) a quasiparticle; jj C — UjCj„ — crvjct a , and the constant Cp 
and the quasiparticle energy Ej are given by 



op = E 



( 6j + ej) - E. 



+ A 2 P /\5 (5) 



and 



1 



2( £ i- £ j)+^' ( 6 ) 



respectively. Here Ap = A(5 J^j ( c ]i c jt) i- s t ne parity-dependent order parameter which has to 
be calculated self-consistently. Ej and Ej are defined as Ej = — [i — \8[u 2 fj S +v 2 (1 — fj a )\ and 

= ^/ Ap + (ej + ej) 2 /4, respectively. For simplicity we neglect the last term of Ej. It does 
not change the qualitative feature of our results since its role is only to increases somewhat the 
effective level spacing near the Fermi level for large 5. fj a = (jj^Jja) is the average occupation 
number of a state j with spin a . It is parity-dependent and given by 

_ U(E j )Z +T f-mZ- 

for P — e (upper sign) and P = o (lower sign), where f±(E) — \j{eP E ± 1) and Z± = 
Iljv(l ^ e~ l3Ej ). Uj,vj are BCS coherence factors: 

•«- 1 -*-i( 1 + Tsf i )- (8) 

The parity-dependent chemical potential /ip is determined by the relation 

(N}= ^log Z P (»)\^ P , (9) 

that holds provided that fip lies halfway between the last filled and first empty level if P = e, 
and on the singly occupied level if P = o ||. 

The mean-field self-consistency condition Ap = X8 Y^ji c ]l c o"t) leads to the relation 

1 = xs E 2f (1_/ ^" /j ' t) ' (10) 

where j c is the cutoff value of j in the summation. At T — with $ = the occupation of 
quasiparticles reduces to fj a — j if the level lies on the chemical potential, namely j — ±jp , for 
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P = o, and zero otherwise. The factor 1/4 is due to the orbital degeneracy of the ring, e° = e?, 
as well as the spin degeneracy. For nonzero external flux, the orbital degeneracy is lifted and 
fjptr = 1/2 (but fj FrT = 0) for odd parity. This nonzero value of quasiparticle occupation gives 
rise to the so called "blocking effect" . That is, the odd parity superconductor has one unpaired 
electron, which prevents pair scattering of other pairs into/out of the singly occupied state 
and reduces the order parameter as compared to that of an even parity superconductor ]To| . 

In the framework of the parity-projected grand canonical description, three different cases 
appear in solving the "gap" equation: (i) Even parity with fully occupied highest level 
(N = Ajp + 2), (ii) even parity with partially occupied highest level (N — 4jp), and (hi) 
odd parity (N — Ajp ±1). It is important to note that the self-consistent equation ( |Io| ) does 
not depend on the flux if (N) is kept unchanged under variation of the flux. The condition of 
constant (N) leads to the flux-dependent chemical potential 

Mf) = M0) + ^ 2 f. (ii) 

Ej is independent of / with this condition, and accordingly Eq.(pO|) is also flux-independent 
at low temperature, fc^T < 5, Ap. 

Fig.|l| shows the pairing parameter as a function of S keeping the electron density constant. 
In solving the equation we chose A = 0.2, close to that of Al p| , and j c = 2jp. The pairing 
parameter which is obtained by solving Eq.(|l0|) depends strongly on the parity as in the 
grain superconductor. In the odd-parity superconductor the order parameter is suppressed 
compared to the one for even parity. However, the behavior of the pairing parameter as a 
function of the level spacing is somewhat different from the grain superconductor with equal 
level spacing || or with randomly distributed levels ||. It is due to the existence of orbital 
degeneracy in the ideal ring geometry. If the highest occupied level is fully filled (case (i)), the 
chemical potential lies halfway between the last filled and the first empty level. In this case 
there exists a critical level spacing S c /Aq ~ O(l) where A e goes to zero as in the mean-field 
solution for the grain superconductivity. If the highest occupied level is partially filled (both 
for even and odd parity) the chemical potential lies on that level. Pair scattering of electrons 
between these orbitally degenerate levels makes it impossible to have a solution with Ap = 0. 
On the other hand, the presence of any weak disorder will lift the orbital degeneracy and give 
a solution with Ap = 0. Note that, in a very small grain with d 3> Ao, quantum fluctuations 
are important and the mean-field treatment becomes invalid |^| . 

Next we discuss the parity-dependent behavior of the persistent current in the "supercon- 
ducting" state. Because we deal with an isolated ring, the current should be calculated in 
the canonical description, while we use the parity-projected grand canonical partition func- 
tion. However, the number fluctuation in the grand canonical treatment is very small at low 
temperature and the free energy with a fixed number N, namely i*jv, can be written as 

F N ~ -~ log Zp + t ipN. (12) 

Note that this relation is exact at T = 0. The persistent current is given by 

I = - c ^¥> (13) 

which gives the following expression by using the mean field Hamiltonian (^) 
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where Idia and I par a are the diamagnetic and the paramagnetic contribution to the current, 
respectively, 

^ - -C S -§N, (15) 

Ipara = — C ^J / ~g^'fj^ • (1^) 

Note that at low temperature (kpT < 8, Ap) the constant term Cp in the Hamiltonian (^) 
does not contribute to the current since it is independent of the flux. 

For even parity, the paramagnetic contribution to the persistent current is absent at fcpT <C 
A e since fj„ = for P = e. Thus, the persistent current for P — e, I e , is equal to Idia and 
obtained from Eq.([ll|) and (|l5|): 

I e = -2/o/ , (17) 

where Iq = evp/(2nR) with vp = being the Fermi velocity and / = ^/^o- (See Fig.^(a).) 
The behavior of I e is essentially equivalent to what is expected in the bulk superconducting 
ring. 

On the other hand, the paramagnetic contribution is important for P = o because of an 
unpaired quasiparticle. Idia is identical for P — o and e. At low temperature (kpT < 8) I pa ra 
can be written as follows 

Ipara — 2Jq (fjpcr - f] F a) ' (18) 

where 

1 e ± ^/* 

f± jF * ^ g (ef/t + e-f/t) " (19) 

t is dimensionless temperature t — kpT/8. At zero temperature I para reduces to sgn[/]io and 
the total current for P = o, 1°, is given by 

7° = / (sgn[/]-2/) . (20) 

As shown in the solid line of FigJ|(b), the current vanishes at half-flux quantum, / = ±1/2, 
while it has maximum value of |/°| at / = 0. This implies that the possible values of fluxoid in 
a ring with odd parity arc half-integers, ±$o/2, ±3$o/2, etc, in contrast to the conventional 
fluxoid quantization at integer multiples of <&o- Note that, the fluxoid is defined in a 
conventional way || 

$' = $ + ^^ s ^l , (21) 

where v s is the supercurrent velocity. 

The parity-dependent behavior of the persistent current found above should be distinguished 
from that of a normal ring of noninteracting electrons with spin 1/2 considered by Loss and 
Goldbart. || In a normal ring with spin 1/2, the persistent current depends on the number 
of particles with modulo 4. Periodicity and oscillation amplitude depend on the number, and 
the current response is diamagnetic only for N = 4jp + 2 (the case of fully occupied top level) . 
On the other hand, a mesoscopic superconducting ring shows only the parity dependence 
that originates from the electron paring. The periodicity and amplitude of oscillation remain 
unchanged as < &o(= hc/2e) and Iq, respectively. 
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The effect of temperature on the persistent current is shown in Fig.^j (b). I para is reduced 
as the temperature is raised. It is clear that the effect of the unpaired quasiparticle disappears 
with increasing T; the persistent current will show a behavior like in Fig.|^(a) , at temperature 
higher than t ~ 1. Actually there is a crossover temperature t* — 0.5 where the current 
response at small flux changes from paramagnetic to diamagnetic. (That is 1 $ (t = t*) = 

o.) 

Let us, finally, discuss the condition under which the parity effect can be observed in 
the persistent current. As noted above, the condition for the existence of the paramagnetic 
contribution to the persistent current is T < T* with k B T* /S = 0.5. Thus, the crossover 
temperature is given by 

k B T* = 2.5(eV)^ , (22) 

where r s is the average distance between electrons. For Al r s — l.I0A= 2.07ao, with ao being 
the Bohr radius. For a ring with R ~ 5 x 10 4 r s ~ 5/xm, T* ~ IK is comparable to the 
transition temperature of bulk Al. If the ring is not strictly ID, the level spacing is reduced 
and, accordingly, the crossover temperature T* is lowered as compared to the ideal ID ring 
with the same radius. Furthermore, T* — > if R — > oo. In other words, the parity effect 
disappears in the thermodynamic limit, and our treatment coincides with the conventional 
BCS description. 

In conclusion, we have investigated the parity-dependent properties of a mesoscopic super- 
conducting ring threaded by a magnetic flux. The properties in the persistent current of an 
even-parity ring are similar to those of a bulk superconductor. On the other hand, a ring with 
odd parity shows unconventional behavior: the fluxoid quantization at half-'mteger multiples 
of the flux quantum. Moreover, there exists a paramagnetic response at temperatures below 
a characteristic temperature of the order of the level spacing. We have also shown that this 
parity effect of the persistent current disappears as the temperature is raised or as the size of 
the ring increases. 

The author wishes to acknowledge valuable discussions with A. Bill. This work has been 
supported by the Visitors Program of the MPI-PKS, and also by the BK21 project of the 
Korean Ministry of Education. 
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Fig. 1. - The pairing parameters as a function of the single particle level spacing for A = 0.2 with 
j c = 2jf- Lines with symbols show the results for three different cases, P = e with N = 4j> (solid 
line with circles), P — o with N = 4jV ± 1 (dotted line with diamonds), and P — e with iV = 4jp + 2 
(dashed line with squares). Ao denotes the bulk gap. 
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Fig. 2. - The persistent current as a function of the dimensionless Aharonov-Bohm flux / = $/<l?o 
with $o = hc/2e (a) for even parity, and (b) for odd parity at t — (solid line), t = 0.1 (dotted line), 
t = 0.3 (dashed line), and t = 0.7 (dot-dashed line). 



